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Abstract. We study the photon blockade phenomenon in a nanocavity containing 
a single four-level quantum emitter. By numerically simulating the second-order 
autocorrelation function of the intra-cavity field with realistic parameters achievable 
in a state-of-the-art photonic-crystal nanocavity, we show that in the strongly coupled 
regime the resulting photon blockade is significantly better than that achievable with a 
two-level emitter. We introduce an intuitive picture of the photon blockade with a four- 
level emitter that explains the performance difference between the two-level and the 
four-level emitter schemes, as well as why - in contrast to a cavity containing a two-level 
atom - signatures of photon blockade appear and should be experimentally observable 
even outside the strong coupling regime when a four-level emitter is used. Finally, 
we show that the emitter-cavity coupling achievable in a nanocavity can overcome the 
non-ideal spacing of optical transitions in realistic four-level emitters that has so far 
prevented experimental realization of this photon blockade scheme. 



PACS numbers: 42.50.Pq, 37.30+i 
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1. Introduction 

In recent years, photon blockade [lj has generated a lot of scientific interest as a 
promising tool for controllably implementing repulsion-equivalent interactions between 
photons, which can be used for quantum simulations [21 El H]. Additionally, photon 
blockade provides a mechanism to coherently generate non-classical light states E] 
for applications such as quantum information processing and quantum cryptography [7J . 
Numerous platforms have been proposed for experimental implementation of photon 
blockade, such as single quantum emitters in cavity quantum electrodynamic (cQED) 
systems [El ED El EE [HI EE2], atoms in Rydberg states [121 [Tlj . and atomic ensembles 
in photonic waveguides [15]. Experimental demonstrations of photon blockade have 
so far been achieved for the most part with cQED platforms in the regime of strong 
coupling. These reports include atomic cQED systems based on a Fabry-Perot cavity 
[9] and solid-state cQED systems based on single quantum dots embedded in a photonic 
crystal nanocavity [H [TB] . Additionally, cQED systems achieving photon blockade in 
the weak coupling regime have been demonstrated with laser cooled atoms coupled 
to a toroidal microcavity [TU] and predicted for quantum dots coupled to either a 
bimodal nano-cavity [H] or to a photonic molecule [T2l ITT] . Unfortunately, platforms 
implementing photon blockade have so far proved to be difficult to scale into more 
complex systems. Photonic crystal cavities containing strongly-coupled self-assembled 
quantum dots have shown a lot of promise in this area [18], but further scalability of 
this platform still remains challenging due to the non-deterministic spatial location and 
spectral properties of the self-assembled quantum dots. On the other hand, cold atoms 
are by nature identical, which makes them excellent candidates for quantum emitters in 
scalable architectures. However, building scalable architecture based on atomic cavity 
QED systems demonstrated so far is non-trivial. Coupling laser-cooled atoms to planar 
photonic-crystal cavities was first suggested in Ref. [19] and has been since explored in 
additional proposals, such as [201 E] • While experimental realization of these proposals 
is challenging, recent developments in design of planar photonic-crystal cavities, such as 
those described in [211 E21 [23] have opened additional prospects of cold atoms interacting 
with the field maximum of a cavity with a mode volume that is ordinarily associated 
exclusively with solid-state cQED. 

In the experimental demonstrations reported so far, photon blockade was achieved 
with a two-level quantum emitter (a neutral atom or a quantum dot) coupled to the 
cavity. However, Imamoglu et al. originally coined the term "photon blockade" in a 
proposal envisioning a four- level atom coupled to a Fabry-Perot cavity pQ . In the strong 
coupling regime, this scheme was predicted to result in a more robust photon blockade 
compared to the conventional approach of a two-level emitter coupled to a single mode 
of a cavity, with quantum interference suggested as a possible cause [21]. In the years 
following the initial proposal, photon blockade based on a four- level quantum emitter 
coupled to a cavity has been analyzed extensively [221 [2S1 [2H [23 [2H1 [2S]- However, these 
theoretical studies examined the blockade under conditions observed in atomic cQED 
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Figure 1. (Color online) (a) Schematics of the described photon blockade system 
consisting of a four level quantum emitter coupled to a cavity. The cavity is driven 
by a weak classical probe with frequency u> p and driving rate £, , while the transitions 
from \g) to |e) and from \s) to |/) are coupled to the cavity field. Levels \s) and |e) 
are coupled by a classical field (drive) that does not couple into the cavity and has 
a frequency u>d and driving rate fl. (b) The bare states of the overall system with 
quantized cavity field and the coupling rates between them, (c) The ground state, 
first, and second manifold of the dressed states of the overall system for the idealized 
atom with A c = A sg = A eg = A/ s = (i.e., probe resonant with the cavity, drive 
resonant with the |e) — > |s) transition, and cavity resonant with \g) — > |e) and \s) — > \f) 
transitions). 

experiments in the regime of strong coupling where the rate of atom-field coupling g far 
exceeds the cavity field decay rate k and, in most cases, only for idealized atoms with 
the transitions \g) — > |e) and \s) — > \ f) (Fig. [l](a)) having the same or nearly the same 
frequency. 

Here, we take a closer look at a four-level quantum emitter coupled to a cavity as 
shown in Fig. [l](a), in the regime of g ~ k with the cavity mode volume V mo d ~ A 3 , where 
A is the wavelength of the probe light. This regime is encountered for self-assembled 
quantum dots embedded in photonic-crystal microcavities, e.g. (301 E], and can also be 
expected for the first generation of experiments coupling laser-cooled atoms to these 
cavities. Note that as an alternative to being trapped at and interacting with the field 
maximum of an air-mode cavity, such as those described in J2TJ [22j [23] , the atom could 
also interact with the evanescent field of a dielectric-mode photonic-crystal cavity while 
being trapped near the surface of such cavity, e.g. using techniques described in [15] or 
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Figure 2. (Color online) (a) Energy splitting of the levels in the first (left) and 
second (right) manifold of the dressed states for an idealized, i.e. resonant, system 
with g eg = gf s = g. Here, Q = g e g/2 (solid lines) and fl = g eg (dashed lines). 
lei 9 — ll 1 9 — 3 x 10 -3 , 7s/ff = 3 x 1CP 4 . (b) Numerical simulation of the second 
order correlation function g^ 2 '(0) of the light transmitted through such ideal system 
for g eg /2ir = 3 GHz and VL — g eg (green curve) compared to that of a two- level system 
with the same parameters (blue curve). The inset plots <r 2 ^(0) on a log scale. For both 
simulations, £/2w = 0.1 GHz and the decay rates of the excited states are 2ir x 10 MHz. 



in j3HE2]. 

As a model system for our numerical simulation, we assume a single quantum 
emitter based on a laser-cooled cesium atom interacting with the field of a solid state 
nanocavity. However, we emphasize that such four-level systems could in principle 
also be implemented in solid state. As the system is driven by two classical fields, we 
monitor the second-order autocorrelation function of the field inside the cavity to study 
the strength of the photon blockade achievable in a cavity with mode volume comparable 
to cubed wavelength of the probe light. 



2. Basic model 



The Hamiltonian describing the coherent dynamics of a system depicted in FigjT|a) 
consisting of a four-level emitter and a cavity is given by [2U 

H = H + Hi n t + Hdrivei (1) 

where, assuming the rotating wave approximation, 



H =u c a Jf a + uj sg \s)(s\ +u eg \e)(e\ +u fg \f)(f\ 

H^t = 9eg{.a ] \g){e\ +a\e)(g\) + g fs ((J\s)(f\ +a\f)(s\) 



Photon blockade with a four-level quantum emitter coupled to a nanocavity 



5 



+ Qe- iuidt \e)(s\ +Q*e iUdt \s)(e\ 

drive 



H HrivP = fe-^'V + S*e itJpt a. 



In this description, a is the annihilation operator for the cavity mode, and the classical 
weak probe field couples to the cavity at a rate S. Another classical field with 
polarization orthogonal to the polarization of the cavity mode provides coupling between 
levels \s) and |e) with Rabi frequency Q without directly injecting photons into the cavity 
mode. The field of the cavity mode couples levels \g) and |e) and levels \s) and |/) with 
rates g eg and gj s , respectively. In the above equations, u c is the frequency of the cavity 
mode and Uij are the frequencies of the atomic transition from level \j) to level \i) in 
Fig. [IJa) and % is set to 1. This analysis is relevant for any single- mode cavity, and does 
not involve any typical characteristics of a nanocavity. By transforming into a rotating 
frame, the original Hamiltonian will change into 

dW 

n = u^Hu + t—u = n + n int + n drive , (2) 

at 

with 

JJ = e _i [ w p( a1 a+\ e ) ( e \) + (u P -Ud)\s) (s\ + (2w p -u d )\f) (f\] 

Ho = - A c a ] a - A sg \s){s\ - A e9 |e)(e| - A /ff |/)(/| 
H int = g eg (a ] \g){e\ +a\e)(g\) + gf s {.a ] \s){f \ + a\f){s\) 

+ n\e)(s\+n*\s)(e\ 
Hdrive &q} -\- S Q, 

where A c = u p - u c , A sg = u p - u sg - u d , A eg = uj p - u eg , A fg = 2u p - u d - u fg . 
In the limit of a very weak probe field {£ — > 0), the couplings between the bare states 
of the lowest manifolds of the overall system are depicted in Fig. [I] (b), where each 
bare state is described by the state of the atom and the number of photons present in 
the cavity mode. The new eigenstates of the system are sketched in Fig. [I] (c) for the 
resonant case (A c = A sg = A eg — Af g = 0). The number of these dressed states states 
in each manifold provides an intuitive explanation of the mechanism behind the photon 
blockade in this system. Specifically, if the probe is resonant with the bare cavity, the 
first photon will be resonant to the \g, 0) — > \D 1: 0) transition, while the second photon 
will be out of resonance with any transition available out of state |-Di,0) [2^1 l25j . If 
we take into account the dissipation of the cavity field to the environment with a decay 
rate k and the spontaneous decay rates between the atomic levels 7 S9 , 7 e9 , 7 es , jf g , and 
jf s , we can describe the dynamics of system with a master equation 

g = -l[n,p] + 2kC[o] + 2 5>y£[<r y ], (3) 

u 

where p is the system's density matrix, ij G {sg, eg, es, fg, fs} and = \ j) C[D] is 
the Lindblad superoperator on operator D used to model the incoherent decays and is 
given by 

C[D] = DpD ] - ^D j Dp ~ \p d]d - ( 4 ) 
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Figure 3. (Color online) (a) A numerical simulation of the second order correlation 
function </ 2 ^(0) of the light transmitted through the resonant system as a function 
of the amplitude of the classical field with respect to g eg for different values of k. 
Here, g eg /2Tr = 3 GHz, £/2tt = 0.1 GHz, gj s = g eg (solid curves), <7/ s = g eg /2 (dashed 
curves), (b) Numerically simulated g^ 2 \0) as a function oigf a /g eg for g eg /2ir = 3 GHz, 
£ /2tt = 0.1 GHz, and Q = g eg /2. (c) Numerically simulated </ 2 )(0) as a function of 
the amplitude of the probe field £ for g eg /2ir = 3 GHz, = g eg . The value of n in 
this simulation is 3g eg (solid blue curve), 2g eg (dashed green curve), and g eg (dash-dot 
red curve). For comparison, the dotted black curve plots gW(0) as a function of £ for 
photon blockade with a two level quantum emitter in a cavity with k — g eg . In all 
three parts of the figure, 7 es = 7 es = 7/ ff = jf s = 2tt x 10 MHz and j sg = 2n x 1 MHz. 



For simplicity, we will ignore additional decay mechanisms and set 7 S = 7 S9 , 7 e = 
leg + les, and 7/ = 7/ 9 + 7/s (Fig. [T] (a)). We numerically simulate the dynamics of 
the system described by equation ^ using the routines provided in the quantum optics 
toolbox [33J with up to six-photon Fock states. 
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3. Fully resonant case 

We start by looking at the idealized case of a fully resonant system (A c = A sg = 
A eg = Af g = 0). This fully resonant case has been explored in detail in previous works 
[2U ESI EHl EH] with parameters corresponding to atomic cQED experiments in Fabry- 
Perot cavities (g eg /2iT = 120 MHz) [34] in the regime of strong coupling (g eg > 5k). Here, 
we revisit the resonant case but perform our numerical simulations with parameters 
that can be expected for a system based on a photonic crystal microcavity, in particular 
Qegjl-n = 3 GHz. This is a fairly conservative value for atom-photon coupling achievable 
in a photonic-crystal nano-cavity, as values of 2ir x 17 GHz have been predicted by Lev 
et al. [20] for closed transitions in cesium. We choose the lower value to budget for 
experimental imperfections, such as the possible use of cavity designs with larger mode 
volume, the use of weaker transitions to implement the four level emitter in a realistic 
atom, and the atom not being localized in the cavity field maximum, perhaps due to the 
atom's thermal motion or because the atom is trapped above the surface of the photonic 
crystal. For the cavity field decay, we focus on the interval of k/2tt « 1 — 10 GHz, 
corresponding to quality factors of Q ~ 3 x 10 4 to 3 x 10 5 . For the chosen value of 
g eg , this range includes the transition point g eg ~ | between the weak and strong atom- 
field coupling in the cavity. Values from this range have been observed in solid state 
cQED experiments in GaAs cavities [5j |30], and values of Q even higher have been 
predicted for air-mode cavities fabricated in other materials [22]. Particular materials 
of interest would include GaP [55] ; SiN, GaN, and Si0 2 [SB], as these - unlike GaAs 
- remain transparent for light at wavelengths corresponding to commonly used optical 
transitions of alkali metals. 

The energy splittings of the dressed states in the first and second manifold are 
plotted as a function of the cavity field decay n in Fig. [2] (a). These plots were 
obtained by numerically evaluating and then plotting the real part of the eigenvalues 
of the Hamiltonian H from equation ^ with complex detunings to include the decay 
mechanisms [27] • For the chosen parameters, the first manifold is symmetrically split 
into three levels, with one of the levels remaining at the energy of the original bare cavity. 
Qualitatively, the behavior of the first manifold is relatively insensitive to the value of 
k and Q. On the other hand, the second manifold splits into three and eventually into 
four distinct energy levels with a decreasing value of k, while the value of k at which 
these splits occur can be affected by changing the ratio between g eg and Q. 

The effects of this level structure on the photon statistics inside the cavity can be 
seen in Fig. [2] (b), where we plot the second order correlation function <?^(0) of the field 
inside the cavity as a function of the gj k ratio with the probe field resonant with the bare 
cavity (A c = 0, green curve). We see that, in a fashion similar to schemes based on a 
two- level emitter coupled to a bimodal cavity pm E] , the photons become anti-bunched 
(<T (0) < 1) aheady for g eg < |, i. e. before strongly coupled regime is reached. This 
is particularly interesting when compared to the energy splittings of the dressed states 
plotted in Fig. [2] (a), since one would not expect a blockade until the energy degeneracy 
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of the bare states is lifted by strong enough emitter-field coupling. For comparison, we 
also plot the g^(0) expected for a two-level atom coupled to a cavity with identical 
parameters and probe detuning optimized for this system (blue curve, A c pa l.bg). The 
anti-bunching resulting from a four-level system is much stronger than that achievable 
in a two- level system, with the difference exceeding two orders of magnitude for g/n > 5 
(Fig. [2^b), inset). This is similar to predictions made based on parameters achievable 
in a Fabry-Perrot cavity [21]. 

In Fig. [3j we numerically simulate the effects of the system's key parameters on 
the photon statistics. For a given value of k, the resulting g^ 2 \0) is fairly robust against 
variations in the magnitude of the coherent field Cl (Fig. [3] (a)). Similarly, p^(0) 
remains nearly constant for gj s > g eg /2 (Fig. [3] (b)). Lastly, the light inside the cavity 
remains anti-bunched for a broad range of values of the probe field amplitude (Fig. [3] 
(c)), before changing into the bunching regime predicted in [29J. 



4. The intuitive picture 



An intuitive picture of why a significantly better photon blockade is expected with a 
four-level emitter can be provided by comparing the transmission of the system for the 
first photon and for the second photon coupled into the cavity. An estimate of these 
transmissions can be obtained by taking the non-Hermitian Hamiltonians characterizing 




Figure 4. Estimated transmission of the first two photons through a cavity with 
a resonant quantum emitter as a function of probe detuning: (a) and (c) for the 
case of a two-level emitter (A es = A c ); (b) and (d) for an ideal four-level emitter 
(A eff = Af s — A c ) with parameters from Fig[2](b). The border of the strong coupling 
regime for the \g) — ¥ |e) transition is marked by the vertical dashed line in each plot. 
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the system's transition between its ground and first, as well as first and second manifolds, 
applying these Hamiltonians to wave-functions describing the system state in these 
manifolds, and then finding a steady-state solution under the condition of a weak probe 
field. 

The non-Hermitian Hamiltonian describing the system's transition between the 
ground state and the first manifold can be written as 

■H {1) = - (zk + A c )a f a- {i lsg + A sg )|s)(s| 

- (i-f eg + A eg )|e)(e| + g eg (a j \g)(e\ +a\e)(g\) 

+ n\e)(s\ +Q*\s)(e\ +£a ] + £*a. (5) 

If we assume the system is in a state described by 

|V (1) > = cM\g,0) + c?\g,l) + cP\e,0) +c^\s,0), (6) 

then the equation we are trying to solve is 

zj t \^ ) )=rt i w ) y (7) 

In the limit of E — > 0, we set c^ ~ 1 and we neglect the manifold-connecting effects of 
the probe, such as the \g,l) — » \g, 2), |e,0) — > |e, 1), \s, 0) — > |s,l) coupling. We define 
the transmission of the first photon though the system as 

r<1> = TT^k " W" 2 ' (8) 

(aTa)o 1^1 

with (a)a)^ being the transmission of the first photon though an empty cavity on 
resonance, and solve equation (7) in steady state (jj<4 = 0) to find the value of c[ . 

In a similar fashion, the non-Hermitian Hamiltonian describing the system's 
transition between the first and second manifold can be written as 

= «W - (z 7/s + A /s )|/)(/| + ^(at| s )(/| + a\f)(s\). (9) 

We assume the system is in a state described by 

\^) = d®\g t 1) + c^\g, 2) + 4 2) |e, 1) + c^\s, 1) + cf|/, 0) (10) 

and we are trying to extract the information about the transmission encountered by 
the second photon (assuming the first photon has already coupled into the system) by 
solving the equation 

Zj t \^)=H^ ] )- (11) 

Again, in the limit of £ — )■ 0, we set ~ 1, neglect the manifold-connecting effects of 
the probe, define the transmission of the second photon as 

T m = J^> „ 2 ^1! k2i (12) 

(aTa)o \ c \ 

with (a^a)^ being the transmission of the first photon though an empty cavity on 
resonance, and solve equation (11) in steady state to find the value of c\ . 
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An analogous approach can be used to estimate the transmission of the first and 
second photon through a cavity with a two-level emitter. The result of that calculation 
is plotted in Fig. [4] (a) and (c) for comparison with the four-level emitter case. We see 
that deep in the strongly coupled regime (k = g eg /3 in Fig. [5] (c)) the best transmission 



of the first photon is achieved for 1 — 1 = 1, while the transmission of the second photon 

9eg 



at this detuning is significantly suppressed (Fig. [4] (a)). Unfortunately, with increasing 
cavity field decay k, the transmission of the second photon increases for |A C | = g eg . One 
can remain in the photon blockade regime by adjusting the frequency of the probe - 
effectively trading off the transmission of the first photon for reduced transmission of the 
second photon - until the transmission peaks of the second photon become too broad 
and adjusting A c will not have a noteworthy effect on the statistics of the transmitted 
light anymore. 

Returning to the ideal four-level emitter, we see that when the probe is resonant 
with the bare cavity (A c = 0, Fig. [4] (d)) the transmission of the first photon is 
maximized, while at the same time there is a minimum in the transmission of the 
second photon (Fig. [4] (b)). This convenient frequency alignment of the maximum 
transmission for the first photon with a local minimum in the transmission of the second 
photon is present even outside the strongly coupled regime and while the energies of 
the dressed states (Fig. [2]^a)) remain degenerate. This minimum can be thought of as 
a generalized case of dipole induced transparency (DIT) [37] and a similar minimum 
in the transmission of the second photon at A c = can be observed in the case of 
a two-level emitter (Fig. [4] (a)). In this case though, this minimum coincides with a 
minimum in the transmission of the first photon (Fig. |4] (c)) and is therefore not useful 
for implementing photon blockade. 

To summarize, the optimal frequency alignment of the maximum transmission for 
the first photon with a local minimum in the transmission of the second photon provides 
an intuitive explanation for the superior photon blockade performance of the four-level 
emitter. 



5. Off-resonant system 

While the fully-resonant system is interesting as a model for understanding the 
mechanisms behind the photon blockade and for estimating the limits of its behavior, 
a practical implementation of this scheme will most likely have to be done with a 
four-level quantum emitter in which the frequencies of the \g) — > |e) and \s) — > \f) 
transitions differ significantly. For example, in the alkali atoms used commonly for 
cQED experiments, such as rubidium or cesium, the frequency difference between these 
two transitions is about two orders of magnitude larger than the atom-field coupling 
rate g achieved in current experimental atomic cQED systems [9j [TOj [38] , which report 
g/27r~ 10 — 100 MHz. The only exceptions are disk cavities reported by Barclay et al. 
[39J, who predict g/2n « 0.9 GHz for atoms displaced 100 nm from the cavity surface. 
However, coupling of atoms to these cavities has not yet been demonstrated to the best 




Figure 5. (Color online) (a) Numerical simulation of g^(0) as a function of cavity 
decay rate k and probe detuning A eg from the \g) — ¥ |e) transition in a four-level 
system based on 133 Cs (ui sg — w/ e = 2ir x 8.941 GHz). Here, the probe is resonant with 
the cavity (A c = 0) and we keep A sg = by setting ujd — u p — tu sg . g eg /2ir = 3 GHz, 
gf s = 1.5g eg , fl — g eg . (b) </ 2 '(0) as a function of the amplitude of the classical 
field il for different values of k. (c) Effects of the probe field amplitude £ on g^(0). 
In both (b) and (c) the value of A eg is chosen to minimize g^(0). In all plots, 
leg = J es = Ifg = Ifs = 2tt x 10 MHz and i sg = 2vr x 1 MHz. 



of our knowledge. 

In this section, we therefore base the four- level quantum emitter on the 133 Cs atom, 
in which u sg — Uf e ~ 2n x 8.941 GHz, to study the photon blockade under conditions 
that are closer a realistic system. Of course this representation of the cesium atom, 
while maintaining the frequency scales, is still highly simplified. To study this system in 
detail one needs to implement a full simulation including the effects of all participating 
levels of the cesium atom in a fashion similar to the work presented by Birnbaum et al. 

Following the approach suggested by Greentree et al. [27J, we keep the probe 
resonant with the cavity, adjust Ud to stay on resonance with the two-photon transition 
\d) ~^ I s ) ( u d = u p — u sg ), and vary the probe detuning A eg from the \g) — > \e) 
transition. In this approach, the probe can couple well into the cavity, while adjusting 
A eg optimizes the strength of the interaction of the cavity photons with the two cavity- 
coupled transitions \g) — >■ |e) and \s) — > \f). Figure [HJ^a) then plots g^(0) as a function 
of cavity decay rate k and probe detuning A eg from the \g) — » |e) transition for a 
probe resonant with the cavity (A c = 0). We see that for each value of k within the 
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Figure 6. Estimated transmission of the first two photons through the system 
numerically simulated in Fig. [5] as a function of detuning A c of the probe from the 
cavity: (a) and (c) for the case of strong coupling regime between the emitter and the 
cavity (k — g eg /3); (b) and (d) for the system being on the verge of strong coupling 
regime (k = 2g eg ). 



range under consideration, there is an interval of A eg for which the light inside the cavity 
becomes anti-bunched (g^ 2 \0) < 1). Similarly, as in the case of the ideal atom, changing 
the value of f2 affects blockade only slightly (Fig. |5](b)) and the light inside the cavity 
remains anti-bunched for a relatively wide range of intensities of the probe light (Fig. 
[5^c)). Note that both of these plots display the minimum value of g^(0) for a given f2 
or S that can be achieved by adjusting A eg while keeping the probe resonant with the 
cavity. 

The effects of adjusting A eg on the photon blockade in an off-resonant system can 
also be observed in the intuitive picture described in the previous section. Figure [6] 
plots the estimated transmission for the first and second photon through the system 
numerically simulated in Fig. [5] (a) for strong coupling regime (n = g eg /3) and outside 
of strong coupling regime (k = 2g eg ). In both cases we see that the transmission of 
the first photon remains highest at A c = as the detuning of the probe A eg from 
the | g) — » |e) transition is changed (Fig. [6] (c), (d)). On the other hand, there is a 
minimum with respect to A eg in the transmission of the second photon when A c = 0. 
This minimum is particularly pronounced in the strong coupling regime (Fig. [6] (a)) and 
remains present, although in a washed-out form, even as the cavity decay is increased 
to k = 2g eg (Fig. [6] (b)). However, terms of the value of A eg at which the minimum 
happens, the agreement with the numerical simulation in Fig. [5]is only qualitative, most 
likely due to the manifold-connecting effects of the probe field £ that are neglected in 
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Figure 7. (Color online) Comparison of g^(0) achievable with a four-level off- 
resonant emitter based on cesium (solid red curve) with a resonant four-level emitter 
from Fig. [2^c) (dashed green), with a comparable two-level emitter (dashed black 
curve), and with a two- level atom coupled to a bimodal cavity (solid blue curve 
and dotted blue horizontal line). The two dotted horizontal grey lines correspond 
to blockade with a four-level off-resonant emitter based on cesium and parameters 
demonstrated in current cQED experiments: g eg /2'K = 120 MHz, k/2it = 40 MHz and 
4 MHz for the upper and lower grey line, respectively. The decay rates of the atomic 
excited states are equal for all plots. 

the intuitive picture. 

The photon blockade achievable with the cesium-based four-level atom is put into 
perspective in Fig. [jj The solid red curve plots the minimum value of g^ 2 \0) from 
Fig. |5](a) as a function of k. We see that the blockade is slightly worse than that 
achievable with the ideal resonant system (green dashed curve), but it still vastly 
outperforms the blockade achievable in an identical cavity with a comparable two level 
atom (black dashed curve). For comparison, the two grey dotted lines mark the expected 
blockade with a four-level cesium-based atom and cavities with mode volume resulting 
in g eg /27r ~ 100 MHz, which roughly represent state of the art cavities used in atomic 
cQED experiments, such as [3IJUIJII2]. The upper line corresponds to k/2tt = 40 MHz 
achieved by Hood et al. [31], while the lower line represents k/27t = 4 MHz, which is 
about five-times smaller than the field decay rate reported by Alton et al. [12] . Finally, 
the solid blue line represents the <?^(0) expected with a two-level emitter coupled to 
two orthogonally polarized modes of a degenerate bimodal nanocavity as proposed by 
Majumdar et al. [TT], with the blue dashed line marking the lower limit of g^(0) in such 
system. We see that this proposal provides photon-blockade performance comparable 
to the four-level atom scheme analyzed here, but its implementation with atoms might 
be difficult due to the selection rules of atomic transitions. 
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6. Outlook 

In conclusion, our results show that, under conditions achievable with currently available 
photonic crystal nanocavities, the originally proposed photon blockade [I] should be 
experimentally observable with realistic four-level atoms and even without achieving 
the strong coupling regime between the atom and the cavity field. The potential for 
scalability of this scheme when implemented with cold atoms and photonic crystal 
nanocavities, as well as the robustness to variations in experimental parameters seen 
in the presented simulations, make the "original" photon blockade a great candidate 
for experimental realization of networks of coupled nonlinear cavities in which the 
interactions between photons can be engineered [2]. 

Additionally, the results presented here might be relevant for exploration of photon 
blockade with four-level solid state emitters coupled to nanocavities, such as nitrogen- 
vacancy centers in diamond [4~3| |4"4"1 14"5] and charged quantum dots in a strong magnetic 
field @B]. 
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